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We show how Geophysics may illustrate and thus improve classical Mechanics lectures concerning 
the study of Coriolis force effects. We are then interested in atmospheric as well as oceanic phenom- 
ena we are familiar with, and are for that reason of pedagogical and practical interest. Our aim is 
to model them in a very simple way to bring out the physical phenomena that are involved. 
LPENSL-TH-06/2000. 



I. INTRODUCTION 

The study of mechanics in non-inertial frames, the 
Coriolis force effects for the rotating Earth being the 
paradigm, are often restricted to the standard examples 
of a deflected projectile and the Foucault pendulum. In 
these two examples, the Coriolis force is only a small 
perturbation. In order to make the Coriolis force effects 
dominant, one must consider phenomena taking place at 
the geophysical scale. 

This is the line followed by the present paper. The 
first section is devoted to the presentation of the basic 
equations. In the second section, we discuss the Physics 
of anticyclones and cyclones: we describe their rotation 
in the geostrophic approximation, and show how cyclones 
(but not anticyclones) may become hurricanes. The third 
section presents a second example of atmospheric circula- 
tion, the Jets Streams, which are stabilized by the Corio- 
lis force. We show also that these strong winds are asso- 
ciated with Rossby waves. The last section presents two 
examples of oceanic circulation: wind-driven circulation 
and Kelvin waves. 



rotating frame TZ', the Navier-Stokes equation takes into 
account the above inertial terms and reads 0: 
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The force / includes the gravitational force and other 
external forces if they exist, p is the density of the fluid 
and p the pressure field. The dependence on M has been 
removed in all the terms for clarity. The centrifugal force 
is conservative. If this is also the case for /, one can 
rewrite the terms Vp, — / and pQ A (i~2 A OM) as Vp', 
where p' is called dynamical pressure. In the rotating 
frame, the hydrostatic equilibrium equation is: Vp' = 
0. The dynamical pressure p' reads, within a constant, 
p' = p + pgz, where g is the Earth gravity field. Recall 
that g includes the centrifugal term, and is thus slightly 
different from the gravitational field, which only takes 
into account the Earth's attraction [pi. 



II. BASIC EQUATIONS 

A. Navier-Stokes equation in a rotating frame 

Let us introduce two frames TZ and TZ' in relative mo- 
tion. The inertial frame TZ is the geocentric one, with 
origin at the center O of the Earth, and whose axes are 
along directions of fixed stars. The frame TZ' is fastened 
to the Earth. It has an angular velocity 17 with respect 
to 1Z, where f2 is the angular velocity of rotation of the 
Earth (ft ~ 7.29 ■ lCT 5 rad • s _1 ). The following relation 
between the accelerations of a point M an(M) in TZ, and 
an' (M) in TZ' may easily be obtained |ij : 



a n (M) = a n ,(M) + 21~2 A v n ,(M) ■ 
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In eq. (2.1), the term 2f2 A Vk'(M) is the Coriolis ac- 
celeration, v-jii(M) is the velocity of M in TZ' , and 
O A (O A OM) is the centrifugal acceleration. In the 



B. Reynolds and Rossby numbers 

The nonlinearity of the Navier-Stokes equation makes 
it difficult to solve in general. It is hence necessary to 
evaluate the relative importance of the different terms 
in order to make further simplifications. This is done 
by introducing the different characteristic scales of the 
flow: L denotes the typical spatial extension, U the ve- 
locity, Q the angular velocity and v the kinematic viscos- 
ity. Two non-dimensional numbers may then be derived 
from these scales. 

(i) The Reynolds number is defined as: 
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It characterizes the relative importance of the momen- 
tum transport in the fluid through advection and vis- 
cous diffusion. For the atmospheric flows studied here, 
typical values are: U ~ 10 m • s — 1 , L ~ 10 km and 
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Thus, the Reynolds number is 
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about 10 10 . A large value of the Reynolds number is 
also obtained for oceanic flows jj). Hence, the Navier- 
Stokes equation reduces, for geophysical flows, to the Eu- 
ler equation: 

— + { V1V ■ V)« w = - - Vp + -f 
at p p 

- 2ft A v n > - ^ A (□ A OM). 

(2.4) 

Moreover, geophysical flows are turbulent (high Reynolds 
number) For the sake of simplicity, we ignore this 
complication in what follows. A simple way of taking 
into account the relevant effects of turb ulenc e will be pre- 
sented in the last section. (See section VA.) 
(ii) The Rossby number is defined as: 
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It compares the advection and the rotation effects. The 
Coriolis force dominates if R <C 1. A geophysical flow, 
characterized by a large spatial extension, may easily be 
influenced by the Earth's rotation, as one typically has 
R ~ 10~ 2 <C 1. On the other hand, an emptying bath- 
tub with U ~ lm • s _1 , and L ~ 10~ 1 m, has R ~ 10 5 . 
Such a flow is more strongly influenced by the advection 
in the fluid, and thus by the initial conditions, than by 
the Earth's rotation. 



III. ATMOSPHERIC EDDIES 

A. Anticyclones and cyclones 

We consider first the situation when the Rossby num- 
ber is negligible. This is the case for anticyclones 
and cyclones since one typically has U ~ 10 m • s -1 , 
L ~ 1000 km, which yields R a ~ 0.1. In the Eulcr 
equation (2.4), we only have to keep the gravity, pres- 
sure and Coriolis terms. This hypothesis constitutes the 
geostrophic approximation. For each point M of the 
Earth, we define a vertical axis (Mz), and a cylindri- 
cal coordinate system (r, 0, z). The vertical component 
of the velocity field is supposed to be zero, which implies 
that the movements of the fluid are locally horizontal, u 
is the radial component of the velocity field and v the 
tangential one. The Earth's angular velocity is writ- 
ten as $7 = On + Oj. where On = f2sinAu z and Oj_ 
is fi's projection on the plane (r, 0); A is the latitude. 
The flow is supposed to be stationary. In this system of 
coordinates, the Euler equation can be rewritten, under 
the geostrophic approximation, as: 
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In these equations, f = 2fl sin A is the Coriolis param- 
eter. In equation ( |3.1c| ), the term 2p(fl± A v) ■ u z i s 
small compared to pg (VLU/g ~ 10 -5 ). Equation (3.1c) 
therefore reduces to the hydrostatic equilibrium equation 
dp/dz = —pg. 

If we consider the case of an eddy in the Northern 
hemisphere and assume that the velocity field is tangen- 
tial (u — 0), then, v < (clockwise rotation) implies 
dp/dr < 0. The pressure is higher at the eddy center 
than outside: it is an anticyclone. A cyclone would cor- 
respond to an anticlockwise rotation. Both situations are 
represented in figure |l|. The rotation senses are opposite 
in the southern hemisphere. 
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FIG. 1. Anticyclone (a) and cyclone (b) in the northern 
hemisphere. 

We end this section with two concluding remarks about 
the presence of the Coriolis force: 

(i) Without this force, an eddy center is always a pres- 
sure minimum. However, in the case of the anticyclone, 
the Coriolis force stabilizes the inverse situation, with the 
eddy center being a pressure maximum. 

(ii) In its vectorial form, the geostrophic equilibrium 
equation reads: Vp' = — 2pfl A v. This implies that 
the pressure p' is constant along a streamline. When the 
usual Bernoulli equation is valid, pressure variations are, 
on the contrary, associated with velocity variations along 
a streamline. 



B. Hurricanes 

Let us consider an eddy (anticyclone or cyclone) whose 
angular velocity and radius are respectively to and R. 
The Rossby number characterizing this eddy can be writ- 
ten as R = U/Lfl = uj/n. Therefore, the geostrophic 
equilibrium corresponds to a small angular velocity of the 
eddy, i.e. oj f2. We shall now consider the case where 
the eddy's angular velocity is not small compared to the 
Earth's rotation. This means that the centrifugal force 
due to the eddy rotation has to be taken into account. 
In this case, the Rossby number is of order unity. In the 
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frame TV , the fluid particle has a uniform circular mo- 
tion. Forces acting on it are the Coriolis force and the 
radial pressure gradient. The equation of motion for a 
fluid particle, located at the eddy's periphery reads, in 
TV: 
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(3.2) 



where ro is the eddy radius. The term —r^u 2 corresponds 
to the centrifugal acceleration of the fluid particle, and 
rofw is the Coriolis term. 

An anticyclone in the northern hemisphere is shown 
in figure |l|a. For such an equilibrium, the Coriolis force 
compensates both pressure and centrifugal forces. If the 
angular velocity of the anticyclone grows, the Coriolis 
force is not sufficient to counterbalance these two forces 
since the centrifugal force grows faster than the Corio- 
lis force with increasing u>. This is not the case for the 
cyclone depicted in the figure |l|b. The pressure and cen- 
trifugal forces may counterbalance each other when the 
rotation of the cyclone becomes faster. This qualitative 
approach shows that there is no limit to the kinetic en- 
ergy of rotation for a cyclone. 



More quantitatively, equation (3.2) can be solved to 
find: 
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where G = dp/dr and Go = prof 2 /4. Figure || gives the 
evolution of an eddy angular velocity as a function of the 
radial pressure gradient. In this figure, the geostrophic 
situation can be found around the origin (small pressure 
gradient and angular velocity). In the northern hemi- 
sphere, the sign of the angular velocity is the same as 
that of the pressure gradient. One can even obtain the 
angu lar velocity of an eddy by developing the expres- 
sion fl3.3j) around zero: u ~ JG/AGo- 
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FIG. 2. Normalized angular velocity as a function of the 
normalized pressure gradient. 

The condition G > —Go, for the existence of the above 
solutions, gives a limit to the angular velocity of an anti- 
cyclone (G < 0). One finds uj max — 2f2sinA. This limit 



does not exist for a cyclone (G > 0). When the angular 
velocity grows, the radial pressure gradient follows this 
evolution and becomes more and more important. This 
explains why hurricanes are always associated with very 
low pressure. 

We note in conclusion that the balance between the 
centrifugal force and the radial pressure gradient is pos- 
sible whatever the sense of rotation. Thus, the existence 
of clockwise hurricanes in the northern hemisphere can- 
not be excluded. However, most of the hurricanes ob- 
served in the northern hemisphere are anticlockwise and 
result from the amplification of earlier tropical cyclones, 
the amplification mechanism being the conversion of the 
latent heat of evaporating tropical warm waters into ro- 
tational kinetic energy. 



IV. JET STREAMS AND ROSSBY WAVES 

The difference in solar heating between the equatorial 
and polar regions drives a convective cell at the plane- 
tary length scale, the Hadley cell, which extends in both 
hemispheres from the Equator up to the sub-tropical lat- 
itudes. The heated equatorial air rises, moves toward 
the poles where it cools, then sinks and comes back to 
the equator. When coming back, the air is deflected 
toward the west by the Coriolis force, generating east- 
erlies at the tropical latitudes which are known as the 
Trade Winds. Conversely, the upper-troposphere tra- 
jectories toward the poles are deflected toward the east. 
Because of the thermal structure of the atmosphere M, 
these upper-level westerlies concentrate in narrow tubes 
of very strong winds up to 80 m ■ s _1 , the Jet Streams. 
The Jet Streams are typically found at altitudes of about 
10 km and at latitudes between 30° and 40°. However, 
their strength and location may depart significantly from 
these mean values, depending on the season, the longi- 
tude, and the day-to-day thermal structure of the atmo- 
sphere at mid latitudes. It can be noted that B. Piccard 
and B. Jones took great advantage of the Jet Streams for 
their recent travel around the world in a balloon. The Jet 
Streams are also useful to the planes flying from America 
to Europe. 

In this section, we propose to show how a zonal wind 
(i.e. along the parallels) may be stabilized by the Coriolis 
force. A mass M of air near the Earth's surface is reduced 
to a point G. Its coordinates are the usual spherical ones 
(R, 6,(p), 9 being the colatitude and R the radius of the 
Earth. The velocity of G can then be explicitly written: 
V H' (G) — R9 ug + RipsmOucp. The quantity Rip sin 9 
is the drift velocity uo of the point G along a parallel. 
We deduce the following expression of the Coriolis force 
moment about the centre of the Earth (point O): 

Mo = 2MR 2 VL9 ug+2MRQ,u cos 9 u v . (4.1) 
The computation of the angular momentum of G about 
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O, in the frame TV , yields: 

L-r,{0) = -(MR 2 )tpsm0u e + {MR 2 )e Uip . 



(4.2) 



The theorem of angular momentum for the point G, 
about O and projected on u v gives: 



v Vu . 
-A = 2 sm A, 
R 



(4.3) 



where A = 7r/2 — 9 is the latitude. This equation is 
linearized for small deviations around a given latitude 
Aq, leading to 



SX + 
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SX = 0, 



(4.4) 



where <5A = A — A . The meridional motion of G re- 
mains bounded, only if uq > 0, which corresponds to a 
drift velocity from west to east. This motion is char- 
acterized by small oscillations around the mean latitude 
Ao with angular frequency ujq = y/ 2f2uo cos Ao / R. These 
oscillations correspond to the stationary case of a Rossby 
wave Q . More generally, Rossby waves in the atmosphere 
are guided by strong westerlies. 



V. OCEANIC CIRCULATION 

Oceanic circulation is, of course, described by the 
same equations as atmospheric circulation. For large 
scale oceanic currents, like e.g. the Gulf stream, the 



geostrophic approximation (see section HI A) is relevant: 
the Coriolis force compensates the horizontal pressure 
gradient, which is related to the slope of the free surface, 
which is not necessarily horizontal [Q. 

We shall be interested here in a slightly different case 
for which the interaction between the wind and the ocean 
gives rise to a current. 



are concerned, the kinematic viscosity is neglected, since 
typically v tU rb/v ~ 10 5 for oceanic flows, and v tU rb/v is 
about 10 7 for atmospheric flows. 

Let us write the Navier-Stokes equation in projection 
on (Oxyz), where (Oxy) is the surface of the globe, (Oz) 
the ascendant vertical, and (u, v, w) are the velocity com- 
ponents: 
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(5.1) 



For a stationary situation with large Rossby number, the 
acceleration terms are negligible: the velocity depends 
then only on space. The horizontal pressure gradient 
terms can also be neglected since the equations have been 
linearized and one can consider the real physical situa- 
tion as the superposition of a geostrophic current (taking 
into account the pressure terms) and a wind-driven cur- 
rent, which will now be described. We consider a solution 
depending on space only through the coordinate z. The 
boundary conditions are the following: the velocity has 
to be finite both as z — > — oo and at the free surface, 
the stress is proportional to dv/dz and parallel to the 
wind flow, assumed to be in the (Oy) direction. One can 
solve eq.(S.l) and find the velocity field (the solution is 
straightforward defining W(;z) = u(z) + iv(z)): 



u(z) = ±V Q cos Q + |) exp (|) 
v{z) = Vbsm ^- 



^ exp V<5 



(5.2) 



where S = ^J2vturb/\f\ is a distance called the Ekman 
depth. Typical values for 5 are 5 ~ 10 — 100 m. "+" 
stands for the northern hemisphere and "— " for the 
southern one. 



A. Wind-driven circulation : Ekman transport 

The wind induces a friction at the ocean surface, trans- 
mitted through turbulence to the deeper layers of the 
sea. There is a supplementary difficulty that we can- 
not ignore here. The flow is not laminar, but essentially 
turbulent. The fluid viscosity is related to molecular ag- 
itation, dissipating the energy of a fluid particle. A dif- 
fusive momentum transport is associated with this phe- 
nomenon. In a turbulent flow, agitation dissipates the 
energy associated with the mean velocity of the current. 
This analogy allowed Prandtl to introduce the notion of 
an eddy viscosity Q. In this approximation, considering 
vjy as the mean flow velocity, the Navier-Stokes equation 
( |2.2D remains unchanged, the eddy viscosity v tU rb being 
added to the kinematic viscosity v. It must be remarked 
that the former is a property of the flow while the latter 
is a property of the fluid. As far as geophysical flows 




FIG. 3. Ekman spiral. The surface is generated by the 
velocity field v(0, 0, z). 
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Close to the surface (z = 0), the current deviates 45°, 
and the direction of the velocity rotates clockwise (anti- 
clockwise) in the northern (southern) hemisphere. The 
amplitude of the velocity decreases exponentially on a 
length scale S, which represents the characteristic depth 
over which the influence of the wind is significant. This 
velocity field, the so-called Ekman spiral, is plotted in 
figure |3j The mean effect of the wind, over a depth 6, is 
the fluid motion in a direction perpendicular to it: this 
effect is called the Ekman transport. 

B. Kelvin waves 



waves is plotted in figure |4|. One can notice that the sur- 
face undulation is trapped in the vicinity of the coast, 
and its spatial extension in the direction of the ocean is 
typically of order L. 




The main difference between atmospheric flows and 
oceanic flows occurs near in the coastline, limiting the 
waters motion. This is the origin of Kelvin waves. If 
one considers the deformation of the free surface of the 
oceans, one can see that gravity acts as a restoring force, 
giving rise to a "gravity wave" ||. When influenced by 
the Earth's rotation, these waves are called "modified 
waves" g. 

Let us consider the following geometry: a south-north 
current, with a coast on its right (east). The coast is 
supposed to be a vertical wall, the water height being de- 
noted ho + h(x,y,t). The Coriolis force usually deflects 
a south-north current toward the east, i.e. toward the 
coast. Hence water gathers close to the coast, and gives 
rise to an west-east horizontal pressure gradient counter- 
balancing the Coriolis force. The equations describing 
the gravity waves are the linearized Euler and continuity 
equations M: 



du 


dh 
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(5.3) 
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dv 

dy 



Taking (Oy) perpendicular to the coast, and considering 
a solution describing the above situation, i.e. v = 0, 
u = £(y) ex P z(<-^ — kx), h = T)(y)expi(wt — kx), one 
obtains: 



u = mo exp 



fy 



exp i(u>t — kx), 




(5.4) 



The dispersion relation is given by u> 2 = ghok 2 , as for 
usual gravity waves. The characteristic length L = 
\fgh~o~l ' f is the Rossby radius of deformation. At a mid- 
latitude A <~ 45°, one finds L ~ 2200 km for ho ~ 5 km, 
while for a shallow sea, i.e. ho ~ 100 m, one rather has 
L ~ 300 km. The surface shape generated by the Kelvin 




FIG. 4. Surface shape generated by a Kelvin wave. 

The Kelvin waves are in fact easily observed, since the 
currents generated by tides are influenced by the Coriolis 
force and give rise to them. As a consequence, the coast is 
always to the right of the flow direction (in the northern 
hemisphere). On the oceanic basin scale, mean move- 
ments are in this case an anticlockwise rotation around a 
point called amphidromic point. This geometry is found 
in many places over the globe, the rotation being clock- 
wise in the southern hemisphere B . 



VI. CONCLUSION 

Coriolis force effects become important as soon as the 
spatial extension of the flow is important (R a ot l/L). 
This is the reason why the Earth's rotation considerably 
influences the atmosphere and oceans dynamics. We have 
presented in this paper several simple examples of geo- 
physical fluid dynamics. We hope it will be helpful for 
Mechanics teachers to illustrate inertial effects with sim- 
ple but physically relevant examples. 
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